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Non-linear  theory  of  alloy  phase  separation  in  open  systems: 
Kinetic  phase  transitions  between  ID  and  2D  structures 

V.  A.  Shchukin  and  A.  N.  Starodubtsev 

Ioffe  Physico-Technical  Institute,  St  Petersburg,  Russia 


Spontaneous  formation  of  nanometer-scale  composition-modulated  structures  is  a common 
phenomenon  for  1II-V  and  II- VI  semiconductor  alloys  [ ].  These  structures  are  mostly 
being  formed  in  open  systems  during  the  epitaxial  growth.  Theoretical  understanding  has 
so  far  been  achieved  of  earlier  stages  of  growth  where  the  linear  stability  analysis  of  the 
homogeneous  alloy  growth  applies  [2-4],  In  semiconductors,  bulk  diffusion  coefficients 
are,  typically,  negligibly  small  compared  to  surface  ones.  Then,  fluctuations  of  composition 
are  created  at  every  surface  atomic  layer  and  are  frozen  once  this  layer  is  overgrown  by 
subsequent  layers.  Frozen  fluctuations  of  composition  from  the  entire  film  affect,  via  long- 
range  strain  fields,  migration  of  surface  adatoms.  This  interaction  can  result  in  a kinetic 
instability  of  a homogeneous  alloy  growth  against  fluctuations  of  composition.  However, 
to  describe  a final  structure  formed  during  the  growth,  a non-linear  theory  is  needed. 

In  our  earlier  paper  [ ] we  have  solved  a non-linear  problem  for  the  final  structure  of  a 
growing  alloy  in  a weak  segregation  regime  close  to  the  onset  of  the  instability,  i.  e.  in  the 
vicinity  of  the  solid  line  at  the  linear  stability  phase  diagram  of  Fig.  1(a).  In  the  diagram  of 


Fig.  1.  Phase  diagram  of  the  alloy  growth.  T is  growth  temperature,  and  v is  growth  velocity, 
(a)  Linear  stability  diagram.  Surface  diffusion  is  isotropic,  (b)  Steady-state  diagram  containing  the 
regions  of  homogeneous  growth,  of  the  growth  of  a 1 D structure,  of  the  growth  of  a 2D  structure, 
and  the  region  where  the  stable  growth  of  both  1 D and  2D  structures  is  possible.  Surface  diffusion 
is  isotropic,  (c)  Linear  stability  diagram.  Surface  diffusion  is  anisotropic.  Directions  of  the  wave 
vectors  of  the  most  unstable  modes  of  composition  fluctuations  are  indicated. 


Fig.  1(a),  at  high  temperatures  T,  the  instability  is  hindered  by  the  mixing  entropy  of  the 
alloy.  At  low  7",  the  surface  diffusion  coefficient  decreases,  D(T)  — Dq  exp(— E^/k^T), 
where  Ea  is  the  activation  energy.  Then  adatoms  are  buried  by  the  incoming  flux  faster  than 
migrate  along  the  surface,  that  hinders  the  instability.  Due  to  elastic  anisotropy  of  zinc- 
blend  semiconductors,  the  instability  occurs  for  the  first  time,  for  the  growth  on  a (001) 
substrate,  for  composition  fluctuations  with  the  wave  vector  in  an  elastically  soft  direction 
[100]  or  [010],  The  solution  of  a non-linear  problem  in  a weak  segregation  regime  | ] shows 
a principal  possibility  of  the  formation  of  either  a ID  structure  modulated  along  [100]  or 
along  [010]  or  a 2D  structure  modulated  in  both  [100]  and  [010]  directions.  However,  a 
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2D  structure  occurs  only  in  a very  narrow  temperature  interval  which  is  hard  to  realize 
experimentally. 

Intending  to  explain  the  formation  of  commonly  observed  2D  modulated  structures,  in 
the  present  paper  we  seek  the  final  state  of  the  alloy  growth  in  a strong  segregation  regime , 
i.e„  in  the  entire  region  of  the  T—v  phase  diagram  where  the  homogeneous  growth  is 
unstable.  We  consider  the  growth  of  an  alloy  A|_cB,C  by  molecular  beam  epitaxy  (MBE) 
on  an  atomically  rough  surface.  The  alloy  is  lattice-matched  on  average  to  the  (001) 
substrate.  The  growth  proceeds  via  deposition  of  atoms  on  the  surface,  surface  migration 
of  atoms  in  a stress-  and  composition-dependent  chemical  potential,  and  incorporation  of 
atoms  into  the  growing  crystal,  desorption  being  neglected.  Let  the  composition  equal 
c(r)  = c + ([> (r ) , the  average  composition  being  c = 1/2,  and  the  surface  profile  equal 
h(x,y ) = vt  + /be,  y),  where  v is  the  average  growth  velocity  controlled  by  the  deposition 
flux.  Coupled  kinetic  equations  describe  a joint  evolution  of  composition  fluctuations  <p(r) 
at  the  advancing  surface  and  of  the  surface  profile  fluctuations  f (jc,  v), 
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Here  F is  the  total  Helmholtz  free  energy,  Df.  is  the  diffusion  coefficient  tensor  related  to 
the  evolution  of  the  surface,  D?  is  the  one  related  to  the  substitutional  diffusion  of  alloy 
components  on  the  surface,  and  a is  the  lattice  parameter. 


Fig.  2.  Coupling  between  composition  modulation  and  surface  profile,  (a)  Two  separate  bulk 
alloys;  (b)  Composition-modulated  structure  with  the  planar  surface;  (c)  Composition-modulated 
structure  with  a non-planar  surface;  (d)  Lattice  parameters  of  bulk  alloys  with  compositions  ci 
and  C2\  (e)  Lattice  parameters  of  coherently  conjugated  composition  domains  in  a stressed  system 
with  the  planar  surface;  (f)  Lattice  parameters  of  coherently  conjugated  domains  in  a partially 
relaxed  system  with  a non-planar  surface;  (g)  Elastic  driving  force  to  phase  separation:  atoms 
A(B ) are  attracted  by  A(B )-rich  domains;  modulation  of  composition  <j>  and  of  the  lattice  parameter 
a;  (h)  modulation  of  the  lattice  parameter  at  the  planar  surface  (dashed  line)  and  at  a non-planar 
surface  (solid  line). 

For  the  lattice-matched  alloy  growth,  the  composition  instability  and  the  morphological 
one  are  not  coupled  in  a linear  regime  [ , ].  A key  property  of  the  non-linear  regime  is  that 
both  instabilities  are  coupled  (Fig.  2).  Let  the  structure  consist  of  alternating  domains  with 
compositions  cj  < c (A-rich  domains)  and  C2  > c (B-rich  domains)  and  the  binary  AC 
have  smaller  lattice  parameter  than  BC.  Then,  for  bulk  alloys  A]_ClBC|C  and  A]_C2BC2C, 
one  has  ao(t’i)  < a = a(c)  < «o(<'2)  (Fig.  2(d)).  In  the  structure  of  Fig.  2(b),  domains 
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are  coherently  conjugated.  Then,  A-rich  domains  are  stretched  with  respect  to  intrinsic 
lattice  parameter  of  bulk  Aj_clBc,C,  i.e.  a(c\)  > ao(ci).  At  the  same  time,  they  are 
compressed  with  respect  to  the  average  lattice  parameter  of  the  structure  a.  B-rich  domains 
are  compressed  with  respect  to  the  lattice  parameter  of  bulk  A|  ,2Br,C,  and  stretched  with 
respect  to  a.  The  surface  consists  of  alternating  domains  under  tensile  and  compressive 
stress  and  is  unstable  against  undulations,  if  the  modulation  wavelength  exceeds  a certain 
critical  value.  Then  the  energy  gain  due  to  elastic  relaxation  exceeds  the  energy  cost  of  a 
non-planar  surface  profile,  and  the  surface  will  consist  of  troughs  over  domain  boundaries 
and  crests  in  the  center  of  each  composition  domain  (Fig.  2(c)).  Due  to  elastic  relaxation, 
lattice  parameters  of  A-rich  and  B-rich  domains  are  shifted  towards  their  intrinsic  values, 
ao(ci)  and  aofe)  (Fig.  2(f)). 

Figures  2(g)  and  2(h)  illustrate  the  reaction  of  surface  undulations  on  phase  separation. 
Atoms  B having  larger  atomic  radius  prefer  to  incorporate  to  surface  regions  which  are 
stretched  with  respect  to  the  average  lattice  parameter  a,  i.e.,  to  B-rich  domains  (see, 
e.g.,  [ ]).  Similarly,  A atoms  prefer  to  incorporate  to  A-rich  domains.  Elastic  driving 
force  to  phase  separation  is  proportional  to  the  actual  variation  of  the  lattice  parameter, 
a ( C2 ) — a(c\).  Since  this  variation  increases  due  to  surface  undulations  (Fig.  2(f)),  the 
elastic  driving  force  to  phase  separation  increases,  too. 

To  obtain  stable  steady-state  solutions  of  kinetic  equations  (1),  we  have  solved  these 
equations  numerically  and  have  checked  the  stability  of  solutions  by  integration  of  kinetic 
equation  over  time  in  the  vicinity  of  the  steady  state  solutions. 

Figure  1 (b)  displays  the  calculated  steady-state  phase  diagram  of  the  alloy  growth  in  case 
of  isotropic  surface  diffusion.  The  modulation  period  is  known  to  increase  with  the  decrease 
of  the  growth  velocity  v [ ].  The  interplay  between  surface  and  elastic  energies  favors 
surface  undulations  at  large  period,  i.e.,  at  low  v.  A non-planar  surface  favors  a 2D  structure 
modulated  in  both  [ 100|  and  [010|  directions  versus  a ID  structure  modulated  in  either  [100] 
or  [010]  direction.  The  effect  is  similar  to  that  for  strained  islands  in  lattice-mismatched 
systems  where  2D  structures  (pyramids)  provide  a more  efficient  elastic  relaxation  than  ID 
structures  (prisms)  [ j.  The  diagram  contains  the  regions  of  the  homogeneous  growth,  of  the 
growth  of  1 D structures,  of  the  growth  of  2D  structures,  and  the  region  where  both  1 D and 
2D  structures  can  grow.  In  the  latter  region,  the  actual  final  structure  of  the  growing  alloy 
depends  on  initial  conditions.  To  obtain  the  structure  for  any  particular  initial  conditions, 
one  needs  to  integrate  kinetic  equations  (1)  over  the  entire  time  of  evolution.  We  do  not 
address  this  problem  here. 

In  Fig.  1 (b)  we  use  typical  material  parameters  for  111-V  semiconductors.  Elastic  moduli 
are  c\\  — l.OxlO12  erg/cm3,  c 12  = C44  = 0.5cn,  the  surface  energy  is  50  meV/A2;  the 
lattice  mismatch  between  pure  AC  and  BC  is  7%;  the  critical  temperature  is  Tc  — 1000  K, 
the  activation  energy  for  surface  diffusion  is  £a  = 1 .5  eV.  Growth  velocities  marked  in 
Fig.  1(b)  are:  v\  — 1.2  A/s,  vi  — 4.0  A/s  and  vj  — 96  A/s.  The  growth  velocity  U3, 
at  which  the  instability  is  suppressed,  is  higher  than  typical  ones  in  MBE.  The  growth 
velocities  ui  and  ih  corresponding  to  the  kinetic  phase  transition  between  the  growth  of  a 
ID  structure  and  the  growth  of  a 2D  structure,  are  just  of  the  order  of  typical  ones  in  MBE. 
Thus,  the  described  kinetic  phase  transition  corresponds  to  realistic  growth  velocities  and 
can  indeed  be  observed. 

A kinetic  phase  transition  between  the  growth  of  a ID  structure  and  the  growth  of  a 2D 
one  has  been  observed  by  Ueda  et  al.  [ ] in  the  MBE  growth  of  InAlAs  alloy  lattice-matched 
to  InP(001)  substrate.  At  the  growth  velocity  v — 3 A/s,  the  growth  temperature  has  been 
varied.  At  440  and  470  °C,  a ID  modulation  in  1 1 10]  direction  is  formed.  The  growth  at 
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higher  temperatures,  500,  530  and  560  °C  reveals  a 2D  modulated  structure  in  [100]  and 
[010]  directions. 

To  address  this  experiment  we  take  into  account  the  anisotropy  of  the  surface  diffusion 
(ZXyjOj  7^  D[iioj)  and  construct  the  linear  stability  phase  diagram  of  Fig.  1(c).  At  high 
temperatures,  close  to  the  onset  of  instability,  the  interval  of  wave  vectors  ky  of  unstable 
fluctuations  is  rather  narrow.  Then,  the  wave  vectors  of  the  most  unstable  fluctuations  kjj0) 
are  parallel  to  elastically  soft  directions  [100]  and  [010],  As  the  temperature  decrease, 
the  interval  of  ky  corresponding  to  unstable  fluctuations  becomes  wider,  the  role  of  elastic 
anisotropy  decreases,  and  the  role  of  diffusion  anisotropy  increases.  The  orientation  of  kj|0> 
deviates  from  [100]  and  [010]  towards  the  direction  of  fast  diffusion  [1 10], 

Finally,  the  direction  of  kj0)  coincides  with  [1 10],  In  this  case,  the  growth  will  pre- 
sumably result  in  a ID  modulated  structure.  Calculated  phase  diagram  of  Fig.  1(c)  is  in 
qualitative  agreement  with  experimental  results.  A detailed  consideration  given  in  Ref.  | ], 
allows  to  reach  a quantitative  agreement  and  to  fit  the  kinetic  phase  transition  tempera- 
ture T\ . 

To  conclude,  the  steady  state  phase  diagram  of  Fig.  1 (b)  applies  to  any  alloy.  We  explain 
a large  variety  of  structures,  e.g.,  the  ID  structure  in  InAlAs  modulated  along  [110]  [ ], 

the  2D  structure  in  InGaAsP  modulated  in  [100]  and  [010]  directions  | ],  and  a kinetic 

phase  transition  in  InAlAs  between  the  growth  of  a 2D  structure  modulated  along  1 1 00]  and 
[010]  and  the  growth  of  a 1 D structure  modulated  along  1 1 10]  | ].  Our  results  demonstrate 
a possibility  to  tune  the  structure  between  ID  and  2D  ones  by  varying  growth  velocity  and 
temperature. 

Acknowledgements 

The  work  was  supported  by  Russian  Foundation  for  Basic  Research  (Grant  98-02-18304), 
and  by  Russian  Federal  Programs  “Solid  State  Nanostructures”  (Project  97-2014)  and 
“Leading  Research  Schools”  (Grant  96.15-96.348).  One  of  the  authors  (V.A.Sh.)  thanks 
the  Alexander  von  Humboldt  Foundation. 

References 

[ 1 ] A.  Zunger  and  S.  Mahajan,  in  Handbook  on  Semiconductors,  ed  T.  S.  Moss,  vol  3,  ed  S.  Ma- 
hajan,  (Elsevier,  Amsterdam)  p.  1399,  1994. 

[2]  V.  G.  Malyshkin  and  V.  A.  Shchukin,  Semiconductors  32,  1062  (1993). 

[3]  J.  E.  Guyer  and  P.  W.  Voorhees,  Phys.  Rev.  B 54,  1 1710  (1996). 

[4]  I.  P.  Ipatova,  et  al.,  Proc.  23rd  Int.  Symp.  on  Compound  Semiconductors,  1996  (American 
Inst.  Phys.  Conf.  Ser.  155,  617);  Phys.  Rev.  B 57,  12968  (1998). 

[5]  V.  A.  Shchukin  and  A.  N.  Starodubtsev.  Proc.  7th  Int.  Symp.  Nanostructures:  Physics  and 
Technology,  St  Petersburg,  p.  497,  1999. 

[6]  F.  Leonard  and  R.  C.  Desai,  Phys.  Rev.  B 55,  9990  ( 1 997 ). 

[7]  V.  A.  Shchukin,  et  al.,  Phys.  Rev.  Lett.  75,  2968  ( 1 995 ). 

[8]  O.  Ueda,  et  al.,  J.  Cryst.  Growth  95,  38  (1989). 

[9]  V.  A.  Shchukin  and  A.  N.  Starodubtsev.  Proc.  Mat.  Res.  Soc.  Symp.  Vol.  583.  Eds.  D.M.  Foll- 
staedt,  A.  Mascarenhas,  and  B.  A.  Joyce,  Pittsburgh.  2000. 

[10]  S.  W.  Jun,  et  al.,Appl.  Phys.  Lett.  68,  3443  (1996). 

[11]  T.  L.  McDevitt,  et  al.,  Phys.  Rev.  B 45,  66 1 4 ( 1 992). 


